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Abstract 

We introduce and solve a semi-classical random walk (RW) model that describes the dynamics 
of spin polarization waves in zinc-blende semiconductor quantum wells. We derive the dispersion 
relations for these waves, including the Rashba, linear and cubic Dresselhaus spin-orbit interactions, 
as well as the effects of an electric field applied parallel to the spin polarization wave vector. In 
agreement with calculations based on quantum kinetic theory [P. Kleinert and V. V. Bryksin, Phys. 
Rev. B 76, 205326 (2007)], the RW approach predicts that spin waves acquire a phase velocity in 
the presence of the field that crosses zero at a nonzero wave vector, go- In addition, we show that 
the spin-wave decay rate is independent of field at qo but increases as (q — qo) 2 for q / qo. These 
predictions can be tested experimentally by suitable transient spin grating experiments. 

PACS numbers: 72.25.-b, 72.10.-d 
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I. INTRODUCTION 



Spin-orbit (SO) coupled two-dimensional electron systems are of great interest, both as 
model systems and as the active component of devices that control electron spin with electric 
fields. 1 Unfortunately, the potential of the SO interaction to control electron spin comes 
with a price - the SO terms in the Hamiltonian break SU(2) spin symmetry. The violation 
of SU(2) means that electron spin polarization is not conserved, decaying instead with a 
characteristic spin memory time t s . The mechanism by which SO coupling leads to spin 
memory loss has been intensively investigated in two-dimensional electron gases (2DEGs) 
in semiconductor quantum wells (QWs), as described in recent reviews.- 1 ^ In GaAs QWs 
and related systems, breaking of inversion symmetry allows SO coupling that is linear in 
the electron wave vector k. 4 ~- The SO terms in the Hamiltonian can be viewed as effective 
magnetic fields that act only on the electron spin, with magnitude and direction that vary 
with k. The loss of spin memory in the effective magnetic field, b(k), takes place through 
the D'yakonov-Perel' (DP) mechanism.-"— In this process the electron spin precesses during 
its ballistic motion between collisions; each time it is scattered b(k) and consequently the 
precession vector, f2(fc), change. The net result is exponential decay of spin polarization at 
a rate approximately equal to Q 2 t, where r is the mean time between collisions. 

There exist two distinct contributions to b(k), the Rashba terrn^ arising from asymmetry 
of the confining potential and the Dresselhaus term^ originating in the intrinsic inversion 
asymmetry of the GaAs crystal structure. A prescription for lengthening spin lifetime in 
QWs of III-V semiconductors by tuning the Rashba coupling strength (a) to equal the linear 
Dresselhaus coupling ((3i) was proposed by Schliemann et at— Recently it was recognized 
that this mechanism amounts to a restoration of SU(2) symmetry even in the presence of 
anisotropic SO interactions.^ The main purpose of this paper is to assess theoretically to 
what extent tuning SO interactions can be expected to increase the distance over which 
electron spin polarization can propagate without decay. 

The potential to extend the spin propagation length despite DP spin memory decay is 
based on the strong correlation between the electron's displacement in space and the rotation 
of its spin on the Bloch sphere. An important step toward a quantitative theory of such 
correlations was made by Burkov et ai— and Mishchenko et a/.,— who derived equations of 
motion that describe the coupling of spin and charge current degrees of freedom in (001) GaAs 
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QWs. Initially only the linear Rashba SO coupling was examined, subsequently Bernevig 
et al. 13 and Stanescu and Galitsk*^ extended the theory to include the linear and cubic 
Dresselhaus terms, respectively. 

The equations of motion can be solved to obtain the normal modes of the coupled system, 
which are waves of mixed electrical current and spin polarization. There exist four such 
modes, reflecting three spin degrees of freedom (S x , S y , and S z ) and the charge density, n. 
For wave vectors, q, parallel to the directions [110] and [110], the four modes decouple into 
two pairs; in one the spin precesses in a plane containing q and the normal direction z, in the 
other the current is coupled to the component of in-plane spin polarization perpendicular to 
Q 

The spin precession mode is the one relevant to spin polarization memory. For example, 
the decay rate of this mode at q = is precisely the DP decay rate, l/r s . In the absence 
of spin-space correlation, the decay rate, 7 ? , of a spin polarization wave would increase 
monotonically with q, i.e., jg = l/r s + D s q 2 , where D s is the spin diffusion coefficient. 
Instead, it was predicted^ that for Rashba SO coupling the minimum decay rate occurs at 
nonzero wave vector, at which point j q is approximately half the DP rate. Bernevig et al— 
showed theoretically that the minimum jg is further reduced when both Rashba and linear 
Dresselhaus interactions are nonzero and vanishes when the strength of the two couplings 
is equal. The resulting "persistent spin helix" (PSH) was shown to be a conserved quantity 
of a newly found SU(2) symmetry that arises when a = (3± and the cubic Dresselhaus term 
(^3) is zero.— However, Stanescu and Galitski 16 showed that perfect SU(2) is broken when 
/?3 7^ 0, leading to large, but not infinite, PSH lifetime. Recently, using the transient spin 
grating technique, Koralek et aZ.— observed the PSH mode experimentally by independently 
tuning the Rashba and linear Dresselhaus couplings. 

The question that arises is whether the PSH effect can be exploited to lengthen the 
distance that a packet of spin polarization can propagate in an applied electric field. In this 
paper we address this question by analyzing the effects of an in-plane electric (E) field on 
the spin-precession modes. We focus on E \\ q, which is the orientation relevant to the drift 
of spin polarization. To predict the spin memory length it is necessary to determine how 
the applied field modifies both the real (9ft) and imaginary (S) parts of the normal-mode 
frequency, oo(q) of spin-polarization modes. The real part is related to the drift velocity 
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whereas the imaginary part is related to the lifetime. The modification of $t{u(q)} is linear 
in E (to lowest order), whereas the affect of E on ^{(^(g)} is quadratic. Kleinert and 
Bryksin— ^ recently have treated this to problem to linear order in E, using quantum kinetic 
theory, and obtained results for 9ft{a;(g)}. 

In this work, we derive and solve equations of motion to quadratic order in E using a 
random walk (RW) approach that is different from previous treatments of this problem. 
The advantages of our approach are physical transparency and mathematical simplicity. We 
construct a semiclassical random walk model that tracks the electron's motion in real space 
and the propagation of its spin on the Bloch sphere. In Sec. [Til we introduce the random walk 
model, derive the equations of motion in the absence of an E field, and solve for the spin- wave 
dispersion relations. We compare the results thus obtained with the earlier quantum kinetic 
theory approaches.— 1 ^ In Sec. IHH we include an in-plane E field, obtaining the equations 
of motion and the dispersion relations to quadratic order. We use the dispersion relations 
to analyze the motion of a spin-polarization packet in the presence of the in-plane field, for 
different regimes of field strength. We illustrate the results by focusing on representative SO 
couplings: linear Dresselhaus coupling only, the SU(2) case where Rashba and Dresselhaus 
terms are equal, and the case of SU(2) broken by a small cubic Dresselhaus term. A brief 
summary is given in Sec. IIVI 



II. RANDOM WALK MODEL 



As mentioned above, as an electron propagates between scattering events, SO coupling 
causes its spin to precess. Thus, as the electron performs an RW in real space, its spin 
performs an RW on the Bloch sphere. We consider a 2D electron gas with both structure 
and bulk inversion asymmetry. The SO Hamiltonian for conduction band electrons in a III-V 
semiconductor QW grown in the [001] direction (taken as ^-direction) is given by, 



H 



so 



fl ■ s, 
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s = hcr/2 is the electron spin, v x and v y are the components of velocity in the [110] and [110] 
directions, a, (3i, and 03 are dimensionless quantities describing the strength of the Rashba, 
linear, and cubic Dresselhaus spin-orbit couplings, respectively, and kp is the Fermi wave 
vector. Spins precess about the effective SO field according to 

ds ^ 

— = Qxs. 3 

dt v ; 

We assume that the impurity potential is short range so that there is no correlation 
between the scattering events. In the absence of the E field, electrons perform an isotropic 
2D random walk with v n (velocity between the nth and (n + l)th scattering events) given 
by Vpt n , where t n = (cos 9, sin 9) is a random two-dimensional unit vector with a uniform 
probability density p n (9) = l/2n. The displacement from nth to (n + l)th step is given by 

r n +\ ~r n = v n r, (4) 

where r is the electron scattering time. In the following we consider Qr, the change in angle 
of the electron's spin between scattering events, as a small parameter. In this case we can 
obtain from Eq. ([3]) the change in the spin direction during the mean-free time as a series 
expansion in Qt, 

As n = s n+1 - s n = n n r x s n + l -a n r x (Q n r x Sn ), (5) 

where we retain terms to second order. 

Let P n {r) be the probability that after n steps of random walk the electron arrives at 
position r and D n (r; s) be the conditional probability that given the electron is at r, its 
spin is s. The joint probability P n (r)D n (r; s) satisfies the following recursion relation: 

P n+1 (r)D n+1 (r; s) = (P n (r - v n r)D n (r - v n r; s - As n )), (6) 

where () denotes average over i n , i.e., (A n ) = J Q 27r A n (9)p n (9)d9. Once P n (r)D n (r; s) is 
determined, the magnetization can be obtained from the following integral on the Bloch 
sphere: 

m n {r)= [ sP n {r)D n {r;s)dZ. (7) 
Js 2 

By substituting Eq. ([6]) into Eq. (J7j), we obtain, 

m n+1 (r) = ( / sP n (r - v n r)D n (r - v n r; s - As n )cE). (8) 
Js 2 



Taylor series expansion on the right hand side of Eq. (jSJ) yields, 

m n+1 (r) =( [ [s + Ct n T x s + \tl n T x (fi n r x s)]{P n (r)D n (r; s 
Js 2 1 



s 2 

- u„r • V[P n (r)D n (r; s)] + \v n T ■ VV[P n (r)D n (r; s)} ■ v n T}dZ). 



(9) 



Again retaining terms to second order, we can write, 

m n+1 = Ji + 1 2 + I 3 , (10) 



where, 



Ji =( / s{P n {r)D n {r; s) - v n r ■ V[P n {r)D n {r; a)] 

Js [ (11) 
+ -v n r ■ VV[P n (r)D n (r; a)] ■ v n r}dT), 

l2 = (f [n n TX8}{P n {r)D n {r;8)-v n T-V[P n {r)D n {r;a)}}dT 1 ), (12) 
Js 2 



and 
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r) = ([ \\il n r x (£l n r x s)]{P n (r)D n (r; s)}dZ). (13) 
Js 2 1 

Upon performing the average over t n , all terms that linear in v n or Q n vanish by symmetry, 
leading to, 

I 1 = m n + Tl op T 2 m n , (14) 



' l \ ny u nx/ \* L nx u ny / q J 1 j 



(15) 



where 



r 2 

^3 = ~— (x(n 2 yn )m nx + y(tt 2 xn )m ny + z(tt 2 n )m nz ) , (16) 

Taking the continuum limit m n — > m(t), (m n+ i — m n ) /r — > dm/dt, and substituting into 
Eq. f JTOj) . we obtain the equation of motion for the magnetization vector. Resolving the 
vector equation into components yields three scalar equations, 

= n °P m * - \{^l) m x - (ttyV x )^^, (18) 
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^^ L = Il o P m y -^(nl}m y + (n x v y }^^, (19) 

= n op m, - l -(n 2 )m z + (n»«.>^ - (20) 

Solving the equations of motion for eigenmodes with wave vector parallel to x yields the 
dispersion relation, 

^ = \ (2(^ 2 ) - m + \{vlW ± yj^ + q^M*. (21) 

This dispersion relation corresponds to modes in which the spin polarization spirals in the x- 
z plane. Note that u)(q) is purely imaginary so that for all wave vectors the spin-polarization 
wave decays exponentially with time. However, the dispersion relation differs from ordinary 
diffusion, where iu oc 1/r + Dq 2 . The difference can be traced to the terms in Eq. (j!5p that 
are proportional to the first derivative of spin density with respect to position - these terms 
are absent in the usual diffusion equation. The coefficients of these additional terms are the 
cross-correlation functions, (Vt x v y ) and (Q y v x ), which shows explicitly that the anomalous 
diffusion is a consequence of the correlation between the electron's motion in real space and 
the propagation of its spin on the Bloch sphere. 

In the SU(2) case (a = j5\ and (3 3 = 0), Eq. ( J2T1) simplifies to, 

«" ± (<7) = \v%r (q ± g ) 2 = D(q ± q ) 2 , (22) 

where D = Vpr/4 and go = ^kpfli. The vanishing decay rate of the cj_ mode at q = 
go indicates the appearance of a conserved quantity - a helical spin-polarization wave or 
"persistent spin helix".— 

The dispersion relations obtained above for the spiral polarization waves are the same as 
those obtained previously, including the cubic Dresselhaus term.— ^ We note, however, that 
while the RW approach accurately describes the spiral coupling of x-z components of spin, 
it does not capture the coupling between charge current and the y component of spin that 
appears in the quantum kinetic formulation. This is because the RW approach does not 
include relaxation to the equilibrium state. In other words, between consecutive scattering 
events the electron's spin precesses about b(k) but has no tendency to spiral in toward it. 
Thus the well-known current-induced spin polarization (CISP) effect 15 is not predicted. To 
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recover CISP requires adding to Eq. Q a phenomenological Gilbert damping term, 

ds 

— = \ G sx(nxs), (23) 
where Xg is the damping parameter. 

III. SPIN HELIX DYNAMICS IN THE PRESENCE OF AN ELECTRIC FIELD 

In this section, we explore how the spin dynamics change in the presence of an E field 
parallel to the wave vector of the spin spiral. To include the effect of E we add a drift term 
to the velocity at each random walk step, 

v n = v F i n + v d x, (24) 

where v d is the drift velocity assumed to be a linear function of E. We assume further 
that the electric field does not change the shape of the impurity potential and therefore the 
scattering probability density is still uniform. 

The drift velocity modifies the precession vector, adding a fixed precession 



Q d = -2yk F 



a + pi 



vp 



v d , (25) 



to Q n at each step of the random walk. Substituting and following the same strategy as 
before, we obtain, 

3th 

I 1 (E) = I 1 -v d T—, (26) 

I 2 (E) = I 2 + fl d rxm, (27) 

h(E) = I 3 , (28) 

where the I\,2^{E) are the quantities Ji,2,3 evaluated in the presence of the electric field. 
The field alters the equations of motion in two ways. First, new terms appear that are 
linear in E. The new term added to I\ converts the time derivative of m to the convective 
derivative, that is the time derivative in a frame moving with the drifting electrons. The 
term added to I 2 indicates that the E field introduces uniform precession about the y axis, 
when viewed in the frame co-moving with v d . The second type of modification is quadratic 
in E; the field increases (fly) by the additive factor Q d and the mean-square velocity (v^.) by 
the factor (v%). 
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Solving for normal modes with wave vector parallel to a;, we obtain 

1 1 //Q2\2 T 2 

iu±(q) = - (2(fi 2 ) - (n 2 )) r + -(vDrq 2 + iv d q ± + {q(Vt y v x )r + ^ d ) 2 . (29) 

To linear order in E, this dispersion relation is the same as that obtained by Kleinert and 
Bryksin.— ^ In the presence of the electric field oj(q) acquires a real part, which describes 
the propagation of spin polarization. Equation ( 1291) also describes the modifications of the 
spin polarization lifetime that appear at second order in E. In the following we discuss the 
spin dynamics that emerge from this dispersion relation for representative SO Hamiltonians. 



A. SU(2) case 

For the case of a = (3i, (3 3 = 0, the dispersion relation simplifies to, 

iu ± (q) =D(1 + 2A 2 ) (q ± q ) 2 + iv d (q ± q ) , (30) 

where A = v^/vp. To distinguish the lifetime and propagation effects we write the dispersion 
relation in the form, 

iu(q) = 7(g) +i<p(q), (31) 

where 7(g) is the decay rate and (p(q) is the rate of phase advance. The real and imaginary 
parts of iu-(q), corresponding to the longer lived of the two modes, are plotted in Fig. [TJ 



As is apparent from Fig. 1(a), the spin polarization lifetime, l/7_(g) remains infinite at 
the PSH wave vector, despite the presence of the electric field. This result is consistent 
with the theoretical prediction that at the SU(2) point the spin helix generation operators 
commute with all perturbation terms that are not explicitly spin dependent.— However, 
the field increases the effective diffusion constant by the factor A 2 so that the decay rate 
for q ^ qo increases rapidly when the drift velocity approaches the thermal velocity of the 
electrons. The spin helix generation operators won't commute with the Hamiltonion if there 
exists a spatial disorder of SO interactions.— 1 ^ 

The rate of phase advance [plotted in Fig. l(b)| vanishes at q — qo, i.e., the PSH is 



stationary, despite the fact that the Fermi sea of electrons is moving by with average velocity 
Vd- Moreover, spin spirals with q < qo will appear to move backward, that is, opposite 
to the direction of electron flow. Although unusual, this property can be understood by 



9 




0.0 0.5 1.0 1.5 2.0 

q/q 



-0.8 

o.o 



0.5 1.0 1.5 

q/q 



2.0 



(a) (b) 

FIG. 1: (Color online) The dispersion relations for (a) the decay rate and (b) the rate of phase 
change of the SO enhanced mode in the SU(2) case, (a) The decay rate 7-(c/) increases with the 
drift velocity (A = Vd/vF) but always vanishes at the resonant wave vector qo. (b) The rate of 
phase change (j)-(q) is proportional to the drift velocity and it crosses zero at the resonant wave 
vector qo. 

considering the spin dynamics in a frame moving with velocity Vd- In this frame E parallel 
to x is perceived as a precession vector fid = —4(3iVdy = —Vdqoy. Therefore in the moving 
frame <f)±{x',t') = ±qx' — Vdqot. Transforming back to the laboratory frame then yields 
0±(g) = v d (q±q ). 

The nature of spin propagation at the SU(2) symmetry point can be made more clear 
if we Fourier transform from the wave vector to spatial domain. If we inject a ^-function 
stripe of z polarized spins at x — 0, the space-time evolution of S z is proportional to the 
propagator, G z (x,t), where 



G z {x,t) oc J dqe iqx (A + e- luJ+t + A-e-*"-*) 



(32) 



where A + and A_ are the weighting factors for the passive and active modes, respectively 
and A + = A_ = 1/2 in the SU(2) case. Upon substituting the dispersion relations uj±(q), 
we obtain, 

f (X ~^ )2 1 (33) 



G z (x,t) oc — cos(gox) exp 



ADt 
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FIG. 2: (Color online) The space-time evolution of S z with a normalized 5-function injection at 
x = 0, t = 0, and drift velocity -u^ = 2Dqo in the SU(2) case. The spin polarization develops into a 
conserved stationary wave with a Gaussian wave packet. 

The spin propagator is the product of a Gaussian envelope function and a static spin wave 
with wave vector go- The envelope function is the one- dimensional diffusion propagator with 
width proportional to y/Dt and drift velocity v&. An illustration of the space-time evolution 
described by this propagator is provided Fig. [2j for a drift velocity = 2Dq . Note that the 
phase of the spin wave modulated by the Gaussian envelope remains stationary as the packet 
drifts and diffuses. This contrasts with the more familiar wave packet, where the modulated 
wave and envelope functions both propagate, albeit with velocities that may differ. 

B. SU(2) broken by cubic Dresselhaus term 

When SU(2) is exact, the integral of the Gaussian envelope function is conserved, even in 
the presence of an E field. However, Stanescu and Galitski 16 have shown theoretically that 
0s, which is nonzero in real systems, breaks SU(2). Koralek et al.— verified experimentally 
that 03 is indeed the factor that limits PSH lifetime in experiments on (001) GaAs quantum 
wells. In this section we calculate the dispersion relation and spin packet time evolution in 
the presence of a small cubic Dresselhaus term. 
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It was shown previously that when (3- 



3 is small, the maximum lifetime occurs when the 



Rashba interaction a = f3\ — [3% (Ref. |16|). We consider a QW with Rashba coupling tuned 
to this value and assume that /?3 This condition is met in QWs in the 2D limit, where 

kpd -C 1 (d is the well width). In this case the dispersion relation in the presence of the 
electric field can be written as 

iu ± {q) = QDklPl + D(q± q ) 2 + iv d (q ± q ) =F iv d Aq, (34) 

where qo = 4kp((3i — j3 3 ) and Aq = 2kpf3 3 . Performing the Fourier transform to obtain the 
space-time evolution of a spin packet, we obtain, 



G z (x, t) oc 7/= e 6jDfc ^3* CO s(g a; - v d Aqt) exp 



(x - v d tf 



(35) 



ADt 

In the presence of the cubic Dresselhaus interaction the integral of the Gaussian envelope is 
no longer conserved. The decay rate can be written in the form, 



3^ 2^3 



2 



illustrating that although the decay rate is nonzero, it is reduced relative to the DP relaxation 
rate by a factor w (f3 3 /(3i) 2 . This ratio is expected theoretically— and has been verified 
experimentally— to be determined by the relation, 

A _ k" F d 2 ( ~ 7 s 

For quite reasonable QW parameters a (3 3 to 0± ratio of 1:100 can be achieved, equivalent to 
a lifetime enhancement relative to the DP spin memory time on the order of 10 4 . 



C. Linear Dresselhaus coupling 

Finally, we consider a fully symmetric well in which only the linear Dresselhaus coupling 
exists. To make comparison with the SU(2) situation, we set the strength of the linear 
Dresselhaus coupling be 2/3 1; so that the resonant wave vector is at q ~ q = Akpfli- The 
dispersion relations j-(q) and <p-(q) obtained by substituting a = f3% = and replacing j3\ 
by 2/3i in Eq. fl29|) are plotted in Fig. [3j Some qualitative features of the dispersion relations 
are similar to the SU(2) case, in that j-(q) has a global minimum and 4>-(q) crosses zero at 
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q ~ go- The most important difference is that the minimum 7_(g) does not reach zero, and 
therefore the spin spiral does decay. In the limit of low electric field, the lifetime of the spin 
spiral is only about a factor of 2 longer than the q = (DP) lifetime. 

The propagation of a spin packet in the linear-Dresselhaus-only case is illustrated in 
Fig. HI using the same initial condition and drift velocity as in SU(2) case. We performed 
numerical integration of Eq. f[3"2"j) to obtain the propagator. As we have seen previously, 
a drifting and diffusing envelope function modulates a spiral spin wave. However, now the 
spiral spin fades very quickly. The contrast between linear Dresselhaus only and SU(2) is 
illustrated in Fig. |5l which is a plot of the integral of the envelope as a function of time. 
After a rapid initial decay, the integral is constant in the SU(2) case, whereas with only 
the linear Dresselhaus interaction the integrated amplitude decays exponentially with rate 
~ Dq 2 . 

Figure M presents another way of visualizing the difference in propagation for the SU(2) 
[Fig. 6(a)| and linear-Dresselhaus-only [Fig. 6(b)| Hamiltonians. The z component of spin 
polarization is shown (with color coded amplitude) as a function of time on the vertical axis 
and position on the horizontal axis. It is clear, from the vertical orientation of the contours 
that the positions of the nodes and antinodes of S z are fixed in space. 



IV. SUMMARY AND CONCLUSION 



We have developed a random walk model to describe the time evolution of electron spin 
in two dimensions in the presence of Rashba and Dresselhaus interactions. From the random 
walk model we derived equations of motion for spin polarization and obtained dispersion 
relations for q parallel to one of the symmetry directions of the Rashba/Dresselhaus Hamil- 
tonian. In Sec. [Til we showed that the dispersion relations for spin-polarization waves that 
spiral in the plane containing the surface normal and the wave vector are identical to those 
obtained from previous analyses.— 1 ^ The random walk approach is instructive in showing, 
in a simple but explicit way, how anomalous spin diffusion and the persistent spin helix arise 
from nonvanishing correlations between the velocity and spin precession vectors. 

In Sec. we obtained dispersion relations for spin-polarization waves that include the 
effects of an electric field parallel to q, to second order in E. The terms linear in E are 
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FIG. 3: (Color online) The dispersion relations for (a) the decay rate and (b) the rate of phase 
change of the SO enhanced mode in the linear-Dresselhaus-only case. The main features resemble 
those in the SU(2) case, both 7_(g) show a minimum and (p-{q) vanishes at qo, but the lifetime is 
finite in this case. 




FIG. 4: (Color online) The space-time evolution of S z in the linear-Dresselhaus-only case with the 
same initial condition and applied E field as in the SU(2) case. The features are similar to those 
in the SU(2) case, except the envelope function decays exponentially. 



14 



(I) 




0.01 



FIG. 5: (Color online) The the absolute value of the spin polarization integrated over position 
as a function of time. In the SU(2) case, (S^l* * is conserved after an initial decay; while in the 
linear-Dresselhaus-only case, (S^l* * decays exponentially. 




(a) (b) 

FIG. 6: (Color online) The space-time images of the spin polarization in the (a) SU(2) and (b) 
linear-Dresselhaus-only cases, respectively. 

equivalent to those obtained from the quantum kinetic approach.—^ To first order in E, 
the field introduces a precession vector in the plane of the 2DEG and perpendicular to E. 
The precession about the y axis gives rise to an unusual behavior in that the spiral with 
wave vector go is stationary in space despite the motion of electrons in the field; waves with 
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q > qo propagate in the same direction as the drifting electrons while those with q < q 
propagate "backward." The terms that are second order in E affect the decay rate of spin 
polarization without changing the velocity. The solutions obtained when these terms are 
included point to the special properties of waves with wave vector q , whose lifetime turns 
out to be unchanged by the field. However, the decay rate of the all other waves increases, 
in proportion to (q — qo) 2 - 

We illustrated these results by considering three representative spin-orbit Hamiltonians: 
SU(2) symmetric or a = (3\ and (3% = 0; SU(2) broken by a small but nonzero P3; and 
linear Dresselhaus coupling only or a = ^ = 0. In order to show the nature of spin 
propagation more clearly, we Fourier transformed the solutions from wave vector to real 
space and obtained the dynamics of spin-polarization packets. In all cases the spin packets 
move at the electron drift velocity. In the SU(2) case the integrated amplitude of the spin 
spiral is conserved, while in the linear-Dresselhaus-only case the amplitude decays with a rate 
~ Dq 2 . When SU(2) is weakly broken by small, but nonzero /?3, the integrated amplitude 
decays at a rate ~ ((3 3 / Pi^Dq 2 ,. 

The conclusions reached by our analysis of the RW model are consistent with a recent 
Monte Carlo study of a specific 2DEG system, a (001) Ini-^Ga^As quantum well with carrier 
density ~ 10 12 cm -2 (Ref. 231 ) . In this study spin polarization dynamics were calculated 
under conditions of steady state injection from a ferromagnetic contact. For a//3i ratios that 
are close to unity, the spin polarization is conserved over several wavelengths of the PSH, 
despite the fact that transport takes place in the diffusive regime. Moreover, the polarization 
is not diminished with increasing electric field. The authors point out that the PSH effect 
can be used to achieve a novel variation of the Datta-Das spin-field-effect transistor (Ref. 



24J ) in which a gate electrode modulates the a to /?i ratio only slightly away from unity. 



This has the effect of varying the wavelength of the PSH without significantly reducing 
its lifetime. Thus small changes in gate voltage can in principle lead to large changes in 
source to drain conductance. Whether such a device can actually be realized depends on 
two factors: fabricating ferromagnetic injectors and analyzers with high figures of merit, and 
demonstrating that the PSH effects that have been observed at temperatures below ~ 100 

n 

K (Ref. 1171 ) can be realized at room temperature. 
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